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We discuss the use of the CP asymmetry parameter (Acp) as a possible observable 
of CP violation in the leptonic sector. In order to do this, we study for a wide range of 
values of L/E the behavior of this asymmetry for the corresponding maximal value of the 
CP violation factor allowed by all the present experimental limits on neutrino oscillations 
in vacuum and the recent Super-Kamiokande atmospheric neutrino result. We work in 
the three neutrino flavor framework. 

PACS numbers: 11.30.Er; 14.60.Pq 



I. INTRODUCTION 



We arc definitively living in very exciting times in neutrino physics. The recent results from Super- 
Kamiokande (SK) indicating evidence for neutrino oscillations in atmospheric showers W , the reports of 
observed neutrino oscillations by the Los Alamos Liquid Scintillator Detector (LSND) M in the — » V e 
and — > v e channels in conjunction with other hints such as the results of solar neutrino experi- 
ments makes it difficult to believe today that all of these facts are not related to neutrino properties 
beyond the standard model. 

In the past neutrino physics has led to the discovery of neutral currents and provided the first indications 
in favor of the standard model of electroweak interactions. It may as well, if neutrino oscillations turn 
out to be confirmed by future experiments, reveal itself as an invaluable tool to cast some light on physics 
beyond the standard model, in particular, on the origin of CP violation. One can hope this can be 
achieved in the study of the neutrino oscillation phenomena in some of the future experiments [Q-|9| . 

The recent KTeV result on Re(e'/e) Jl0| finally establishes direct CP violation in the Kaon system, 
rules out once and for all pure superweak theory and supports the notion of a nonzero phase in the CKM 
matrix. This makes it even more interesting to check if a similar effect also happens in the lepton sector. 

As it is well known Jll]-|l4|] CP violation in neutrino oscillations can, in principle, be observed in 
neutrino experiments by looking at the differences of the transition probabilities between CP-conjugate 
channels, AP = P(V a — > Vp) — P(v a — ► vp). It has been pointed out by many authors that it may be, in 
practice, very difficult to get a reliable measurement of AP due to possible Earth matter effects in long 
baseline neutrino experiments p^ , |l6| . However we believe it is worthwhile to try to evaluate the maximal 
size of CP violation effect in vacuum since the vacuum oscillation experiments will be the ultimate proof 
of electroweakly induced neutrino flavor conversion. Both solar and atmospheric neutrino data can, in 
fact, be explained by alternative mechanisms invoking matter phenomena []1 T] — pL 9|] - 
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Here we will use the asymmetry parameter, Acp — AP/[P(v a — > vp) + P(y a — > vp)], suggested 
by Cabibbo Q, as an alternative to measure CP violation in the leptonic sector. We will investigate 
the possible values of this parameter and the corresponding maximal values of AP allowed by present 
experimental data for different L/E situations, for a particular choice of neutrino mass squared differences. 
We will not make any assumption about the elements of the mixing matrix. 



II. SCALE OF MASSES AND CP VIOLATION PARAMETERS 



This work will be developed in the three flavor neutrino scheme and for this reason only two mass scale 
indications can be taken to be right (20). We will fix them to be: 

Am 2 21 » 3.0 x f(T 3 eV 2 , Am\ 2 w 0.4 eV 2 , (1) 

which are taken within the allowed regions given by the atmospheric |lj and terrestrial neutrino experi- 
ments ^J2l],|22| respectively. We will also take into account the probability constraints coming from this 
two types of experiments. We admit here that the solar neutrino problem may be understood invoking 
other types of mechanisms |l7],p3| . 

We can get the analytical expressions for AP and Aqp using the usual form of the CKM matrix 
parameterization: 



U = 



C12C13 S12C13 Sl3e~ 4<5 

-S13C23 - Ci 2 S 2 3Sl3e lS C12C23 - Si 2 S 2 3Sl3e lS S23C13 
Sl2S23-Cl2C 2 3Sl3e l5 -Cl2S23-Sl2C 2 3Sl3e I<5 C 23 Ci 3 



(2) 



where c and s denote the cosine and the sine of the respective arguments. 
Thus AP in vacuum can be written as: 

AP(a, 0) = P{V a -> Vp) - P(v a -» vp) = 4J C p(sin A12 + sin A 23 + sin A31), (3) 

with a, f3 = e, /j,,t and 



Am 



leV / \ km / V E J 
where Amf- = m| — m? and the well known Jarlskog invariant p4| 

J CP = C 13 Si3Ci2Sl2C23S23Sin5. (5) 

We can see from Eqs. (^-||) that in vacuum AP(/j,, e) = AP(/i, r) = AP(e, r) so they are all simply 
referred to as AP. 

On the other hand A(a, 0)cpi which depends on the specific channel (a, (5), is given by : 

A( m _ P{v a -> Vp) - P{v a -> vp) , . 

A{a,l3) C p - \ , or \' ^ b 

P(v a -> vp) + P(y a -> vp) 

In practice, for a real experiment, one has to average the probabilities in Eqs. (^) and (g) over the 
corresponding experimental conditions, i.e. neutrino energy spectrum, distance L, efficiencies etc. 



III. EXPERIMENTAL CONSTRAINTS 



We have studied the experimental constraints imposed on the neutrino oscillation probabilities by sev- 
eral experimental results. In order to obtain realistic constraints we have, for each experiment, averaged 
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the corresponding probability over its distributions of L and neutrino energy. This will be denoted by 
(• • •) in what follows. 

In the terrestrial case we have taken the results from two short baseline accelerator experiments 
LSND § and E776 || for the channels Vu — ► v e and — > v e respectively. To analyze LSND we 
have used the simple model quoted in Ref. p5| . In the case of E776 the neutrino beam energy spectrum 
was taken from Ref. |p2fl , 

LSND: 

(P(i^ V e )) = 0.31 ± 0.10 ± 0.05%; L ~ 30 m, E ~ 36 - 60 MeV, (7) 



E776: 



(P(^ -> v e )) < 1.5 x 10~ 3 at 90% C.L.; L/E ~ 1. (8) 



For the oscillation channel — > i/ T we have used the data from CHORUS [£6| . The CHORUS relevant 
information is available at their Web site [^7|. This short baseline accelerator experiment gives 

(P(v„ -> v T )) < 6 x 10~ 4 at 90% C.L.; L/E ~ 0.02. (9) 

We have obtained from CHOOZ |^8| , a long baseline reactor experiment searching for the disappearance 
of the 77 e , the limit : 

1 - (P(v e ->F e )) < 10" 1 at 90%C.L.; L/P- 300, (10) 
and from the Bugey |pl[| reactor experiment 

1 - (P(V e -> F e )) < 10~ 2 at 90%C.L.; L= 15, 40, 95 m, B~ 1-6 MeV. (11) 

We have read out the positron energy spectrum for CHOOZ from Ref. [^8) and for Bugey from Ref. [2lf| 
and used the relation E v = E e + +1.8 MeV to extract the neutrino energy spectrum. 

We show in Fig. |l| our own reproduction of the exclusion (allowed in the case of LSND) regions of 
these experiments. We observe a reasonable agreement with the original experimental plots that can be 
found in Refs. P, pTp^j2^ , |28t . This is a demonstration that Eqs. correctly represent the constraints 

coming from these reactor and accelerator neutrino experiments. 

In the case of atmospheric neutrinos, we only have used the latest SK Q result since it is the most 
precise one. We have performed our analysis of SK atmospheric data in the following way. We have 
separately used the muon neutrino events (/i-like events) ratio 



R,{L/E) = 7 ({P^{L/E)) + ^(P eii (L/E))) 
\ Mo / 

and the electron neutrino events (e-like events) ratio 



(12) 



R e (L/E) = 7 ((P ee (L/E)) + ^(P^L/E))) , (13) 

where fio = /-Iq(LJE) and eo = eo(L/E) are the distributions of Monte Carlo //-like and e-like events 
taken from Ref. [ |30[ and 7 is a constant to take into account the overall normalization. We have taken 
the value of this constant to be 1.16 in accordance with Ref. |^|. We have introduced in Eqs. ( fl2"| ) and 
( p^ ) the shorthand P a p — P(y a — * ^),ct,/3 = e,/i. The probabilities (P) have been smeared by the 
resolution function given in Ref. |u| at each L/E value. 

On the other hand, from Fig. 4 of Ref. § we can read eight different values of L/E and their corre- 
sponding /i-like and e-like ratios. Using this we have calculated the simple average of each ratio obtaining: 

=0.76 ±0.08; R e = 1.19 ±0.13, (14) 
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defining two allowed bands, one for the e-like events and the other for the /x-like events ratio. These are 
our selection criteria coming from SK data which were used in the following way: for a given set of values 
of the mass squared differences, the mixing angles and the phase we have computed the right hand side of 
Eq. ([12]) and E q. (fl3| ) taking an average over the eight L/E bins and verifying if they satisfy the conditions 
given by Eqs. (|14|). It is important to stress that working with /x-like and e-like event ratios separately 
has the advantage that we are able to constraint independently P(f M — > v e ) and P(y e — > v^). Indeed 
doing this we observe that these probabilities are always below a few %, this is much more stringent that 
the SK analysis shown in Ref. |3^j3^ |. 

It must be noted that this simple SK analysis is not totally rigorous, in fact by using the L/E distri- 
bution from Ref. |Q we only work with a sub-sample of the SK data (fully-contained events). Moreover 
there are quite large uncertainties in the determination of this distribution due to the fact that one relies 
on the observed final lepton to infer the neutrino physical quantities, these uncertainties are discussed in 
Ref. Nevertheless we consider our approach good enough for the goals of this paper. 

since the Eqs. (00) are in fact independent of the number of neutrino generations we are now free to 
use these probability limits in the three neutrino flavor framework. 



IV. LIMITS OF J cp, AP AND A C p 



In order to determine the maximum permitted values of AP, A(jj,, e)cp and A((i, t)qp we proceed in 
the following way. We choose randomly different values of the mixing parameters s| 2 ,s|3 and s 2 3 taken 
in the interval [0, 1]. For each drawn set of these parameters we evaluate the values corresponding to Eqs. 
(7[-0[). We check if they simultaneously pass all the experimental constraints given in Eqs. (0-|ll|) and Eq. 
(14|)T It really means that we look for a common allowed region among all experiments. In the positive 
case we compute the CP violation factor Jqp for fixed values of the phase angle 8. We select among them 
the maximal value of Jcp, Jcp^i f° r each 6 and calculate the corresponding values of AP max , A(/j,, e)cp 
and A(h,t)cp for L/E varying from 40 to 250. 

In Tables ||| we show J™|, x for Am^ w 3.0 x 10~ 3 eV 2 and Am§ 2 « 0.27, 2.0 eV 2 . We have picked 
these values so one can have an idea of the order of the variations that J™p x will have if one varies Am| 2 
inside the LSND allowed region. In general we see that the values of J™p x in the squared mass difference 
interval considered in this paper are more or less stable and of 0(1O~ 3 ) in agreement with Ref. |M. 

In Fig. ^ we show AP max as a function of L/E for sin 6 = 1. As one could guess the oscillating behavior 
of this curve is dictated by the sum of sines in Eq. (J3J) - We are not showing here other curves for different 
values of sin 5 since they have exactly the same form as these ones, only the oscillation amplitude will 
change in accordance to the corresponding J™p x . We observe that the maximal value of AP max , which 
increase with L/E, are of the O(10 -2 ). This is in accordance with the estimations given in Refs. | ^6| , p4[ . 

In Figs. HH we show the behavior of A(/j,, e)cp and A(fi, t)cp as a function of L/E for sin 5=1. We can 
observe that A([i,e)cp grows and A(fi,r)cp decreases as a function of L/E. Also the maximal values 
for A(fi,e)cp are of the 0(1) while for A(fi,r)cp they are of the O(10 _1 ). This result was expected 
since we know that the only difference between A(fi,e)cp and A(fi,T)cp, in vacuum, comes from the 
denominator of Eq. (||) and P(v^ ~ > v e ) is currently very much suppressed by data while SK data seems 
to support P(v^ — > v T ) oscillations that can reach the order of 10 _1 for neutrinos coming from bellow 
the horizon. 

As we have already mentioned at the end of Sec. || one has to be careful in interpreting our results 
in relation to future experiments and remember to take into account the average over the corresponding 
L/E distribution. Because we did not want to make any ad-hoc hypothesis on L/E distribution we 
have deliberately chosen not to do any smearing on L/E. We believe that presenting our results in this 
unfolded way make them more useful and ready to be applied to any real experimental situation. 

To illustrate the implications of our results in future experiments we have put an arrow in Fig. 0^ at 
the corresponding mean value of L/E for MINOS §, K2K § and one of the possible configurations of 
a long baseline neutrino experiment using neutrinos from a muon collider p5| . Based on our results we 
also have computed an estimation, averaging over the expected energy spectrum, of the maximal values 
that can be investigated at MINOS: (A(^,e) C p) = 0.33 , (A(fj,,r)cp) = 0.02, (AP) = 0.0022; and K2K: 
(A(fi,e) C p) = 0.19 , (A(h,t) C p) = 0.0017, (AP) = 0.0014. All these values are, as expected, lower 
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than the ones shown by the corresponding arrows in Fig. Q-^ since in the L/E scope of MINOS/K2K the 
dominant scale is Ato 2 ,-^ the contribution from Am§ 2 being averaged out. 

V. CONCLUSIONS 

We have found the maximal allowed values for AP max , A(/z, e)cp and A(fx, r)cp as a function of L/E 
and sin 5. This was done in the three neutrino flavor framework using the most stringent constraints from 
recent neutrino data and admitting the two mass squared differences to be Ato 2 ^ ~ 3.0 x 10 -3 eV 2 and 
Ato§ 2 ks 0.4 eV 2 . In fact this is in a way complementary to Ref. 6q ]. 

It is important to remark that we have adopted here a different approach from the authors of Refs. 
|l^,|34|,|36| since we do not make any assumptions about the mixing parameters other than the two squared 
mass difference scales and work with the probability expression without any approximation. Besides we 
have explicitly used the experimental resolution functions in our calculations. 

We have seen that in general the values of Aqp are much more sizable than the corresponding AP max . 
In fact, admitting the mass hierarchy considered in this paper, we have shown that the present neutrino 
data tremendously suppress the maximal values AP that can be investigated at the next generation of 
neutrino experiments. In addition since in the asymmetry the systematic errors cancel out, even if the 
absolute flux of the neutrino beam is determined with an accuracy of 10 % it may be possible to measure 
CP violation at 1 % level. This is particularly interesting since in long baseline experiment there are 
expected matter effects that fake genuine CP violation . In a forthcoming paper |5?J we will discuss 
the implications that the inclusion of matter effects will have on our present limits as well as how these 
limits change if one varies the squared mass difference scales. 
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4 x 


4 x Jgf7sin<5 


sin 8 


0.0011 


0.0044 


0.2588 


0.0024 


0.0048 


0.5000 


0.0032 


0.0045 


0.7071 


0.0046 


0.0053 


0.8660 


0.0052 


0.0054 


0.9659 


0.0052 


0.0052 


1.0000 



TABLE I. Maximal values of the Jarlskog factor obtained for different values of sin 5 with 
Amir ~ 3.0 x 10~ 3 eV 2 and Am| 2 ~ 0.27 eV 2 . 



G 



4 x Jgf? 


4 x Jg??/sm6 


sin 8 


0.0010 


0.0041 


0.2588 


0.0021 


0.0042 


0.5000 


0.0031 


0.0044 


0.7071 


0.0038 


0.0043 


0.8660 


0.0042 


0.0044 


0.9659 


0.0043 


0.0043 


1.0000 



TABLE II. Maximal values of the Jarlskog factor obtained for different values of sin S with 
Amli « 3.0 x 10~ 3 eV 2 and Am'L « 2.0 eV 2 . 
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FIG. 1. Our reproduction of the exclusion/allowed regions for different experiments in two generations: (a) 
E776, (b) LSND, (c) Bugcy, (d) CHOOZ and (e) CHORUS. The exclusion regions at 90% C.L. are the ones to 
the right of the curves for all experiments except LSND. In the case of LSND the allowed region is the region 
between the two curves. 
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